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1. Introduction 



Quantum chromodynamics (QCD) is the theory describing the strong interactions. Accord- 
ing to QCD, as the temperature increases, hadronic matter undergoes a transition to quark gluon 
plasma. The function, which describes the equilibrium in this system for different temperatures, is 
the equation of state. For the equation of state, practically one has to measure only the pressure, 
which then will unambiguously determine the energy density as well. A difficulty emerging about 
the equation of state is, that perturbation theory does not seem to fit lattice results even at higher 
temperatures (smaller couplings). While available lattice results both for pure gauge theory [[I], 
and for full QCD j3|, |], ||] (and see also ^) end at around 5 • Tq, standard perturbation theory 
converges only at extremely high temperatures. To create a link between these two methods, we 
present lattice results on the pressure at temperatures, which were previously unreachable. It be- 
comes possible to compare our data with perturbation theory formulae. Our results are obtained 
using two different approaches: we present a new way to renormalize the pressure (for this we have 
results for N t = 4), and a direct method to measure the pressure (results for N t = 4, 6 and 8). 



2. Renormalization of the pressure 

For illustration, we present here the technique and the results for pure SU(3) gauge theory. 
The extension to full QCD is straightforward. 

In order to obtain results which can be extrapolated to the continuum, one has to renormalize 
the pressure. This is usually done using the standard integral method [Q] : 

p^{T)=p{T)-p(Q) = |^((Pl) r -(Pl) ) (2.1) 

This way one has to carry out simulations on finite and zero temperature lattices. However, exactly 
zero temperature cannot be realized. Together with the fact, that divergences that are removed by 
renormalization are independent of the temperature, rises the question: why not use finite T lattices 
for renormalization? In order to do so, let us introduce the following quantity: 

p(T)=p(T)-p{T/2) = jd$ ((PI)*, - (Pl) 2 . Ni ) (2.2) 

Here, instead of using 1/2 for the subtraction temperature, one can use any factor, which is smaller 
than 1. Now, we may build up p reB as a sum of differences, so it can be expressed with p as: 

p m (T) = p(T) - p(T /2) + p(T /2) - p(T /4) + ...= p(T) + p(T/2) + ... (2.3) 

In fact, we usually measure the dimensionless pressure, which can be obtained by including the T 
factors: 



Pren _ P_ 
rp4 jt4 



T 16 r 4 



T/2 256 



Due to the increasing powers of 1/16 in the forthcoming terms, one practically needs only a few of 
them. Applying this scheme, we may reach arbitrary high temperatures using lattices with only N t 
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and 2 • N t temporal extent. It is worth mentioning, that a similar formula can be constructed for the 
case of the normalized interaction measure / = (e — 3 • p)/T 4 : 



where 



I im (T)=i(T) + ^I(T/2) + . 



I(T) = N? ■ dp/dloga ■ ((Pl>,v, - (Pl) 2 .iv ( ; 



(2.5) 



The above method can also be easily generalized for the case of dynamical fermions. 

We can demonstrate the applicability of this technique by using only the finite temperature 
N t = 6 and N t = 8 data of |[T]] to reproduce their results obtained by the standard method (which 
uses the more expensive T = data). Here we will use the ratio x = 6/8 instead of the previously 



shown 1 /2. We keep only four terms in the sum (2.3), since the forthcoming terms are practically 
negligible: 



Pren 



p{T) +x A ■ p(xT) +x 8 • p{x l T) +x lz ■ p{x 5 T 



12 



(2.6) 



On figure |]. we show the results obtained by the standard renormalization procedure, and also 
those by our method. 



— i i i i i | 

■ Boyd et al. N,=6 

— final- sum of p' 

— p' 




Figure 1: The dimensionless pressure as a function of the temperature. By summing up the intermediate 
terms (indicated by solid blue lines, denoted by p 1 ) we arrive at the total renormalized pressure (solid red 
line), which agrees completely with results obtained by the usual integral method. 



3. Setting the scale 

Our new method is not only capable of reconstructing results of the standard procedure, but 
it also has a clear advantage over that, since T = lattices are not needed for the subtraction, 
only for scale determination. So we need large statistics only for T > simulations, which are 
less demanding in terms of computer resources (memory, CPU time). If we wanted to reach very 
high temperatures, we still needed large, and therefore expensive zero temperature lattices to set the 
scale. Alternatively, in the asymptotic scaling region, one might determine the scale using improved 
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perturbation theory. To verify that we are indeed in the asymptotic scaling region, we compare 
results of the Sommer parameter from [JsJ] , and a 3-loop improved perturbation theory formula [§] 
inserted into the expression of the lattice spacing a = A~ l (a — > 0)/i mp pr(j8). This comparison on 
figure^] shows that the last few lattice simulation points are in the asymptotic scaling regime, which 
means that an extrapolation with the perturbative formula might be trusted. 




Figure 2: Asymptotic scaling seems to be realized, as the improved perturbation theory formula fits the 
lattice results of the Sommer parameter for larger values of /3 . 

4. Results I. 

For our simulations we used tree level improved Symanzik action, and an overrelaxation- 
heatbath algorithm. The simulations were carried out on lattices with temporal extent N t = 4; we 
performed the subtraction on N t = 8 lattices. In order to account for decreasing screening masses, 
we used a rather large aspect ratio of N s /N t = 8. On figure ||[ we show our results on the pressure. 




T/T c 



Figure 3: We show the two, non-negligible terms (denoted by dotted lines) of the normalized pressure, and 
their sum (solid line). P$b denotes the pressure of the non-interacting gluon gas. 

It is worth mentioning here, that independently of the renormalization procedure, a problem 
emerges within the integral method framework. Since strictly speaking, the pressure is only exactly 
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zero at T = 0, in principle we would have to cany out the integration stalling from zero temperature. 
Due to the uncertainty of setting the lower point of the integral, and also due to larger statistical 
fluctuations in the low temperature regime, we needed about 6 times more statistics for the region 
below 7c, as we did for T = (1 . . . 100) • Tq. 



5. Direct approach 



In the following, we will present a method, which is suitable for measuring the pressure di- 
rectly, instead of integrating the difference of two, separately measured plaquette variables. This 
way we can get rid of the above mentioned lower point-related uncertainty in the integral: we can 
provide a reference point in the pressure. 



Let us consider one single term of (2.3), namely p as defined by eq. (2.2): 



1 



■logZ(N t ] 
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■logZ(2AT f ) 
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(5.1) 



y N,Nf 6 v " 2N t N* b v ' 2N t Nf ™° \Z(2N t ) 
where we expressed the pressure as logZ. Now schematically we can draw the ratio of the two par- 
tition functions like the following (let S 2 b be the action for the boundary condition in the numerator, 
and Sib f° r the one in the denominator): 

N t -2 N,-l 



Z 2 (N t ) 
Z(2N t ) 



Now let us take an interpolating partition function Z(a) = J £>Uexp(—(a ■ S2b + (1 — a) • Sib)), 
which could be depicted like: 

a r (l-o) 

Z(a) = 



Using Z(a), one obtains 



log 



lot 



Z(V 



(l-«) 



1 , <flogZ(a) 
da — 
o da 



da{Sy h -S 2 b) 



(5.2) 



\Z(2N t ) J " e \Z(0) j jo aa jo 

So we can calculate the pressure itself at any given temperature, without carrying out simulations 
at lower temperatures, and then performing an integral. However, the cancellation of the plaquette 
variables, which makes the former method a difficult task, also emerges in this new method. On 
figure we plot the integrand of (5^) as a function of a to show why it is hard to calculate the 
integral with a given precision. 

Still, it is worth using this new method, since our integrand here is proportional to N~, so it is 
expected to scale with Nf instead of ./V f ~ 4 - This means that we gain a factor of with respect to 
the standard method. In fact, we have to perform the integral in a for each j8 we need. However, 
there is no need to worry about the lower part of the temperature interval anymore, since using this 
method we are able to set the integration constant in the pressure. 
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Figure 4: Cancellation effect in calculating the pressure by our new method. These results were obtained 
on N t — 4 lattices at = 50 (which roughly corresponds to the Planck temperature). 



6. Results II. 

The method discussed above gives us the possibility to measure the pressure at very high 
temperatures. This was carried out using lattices with temporal extension N t = 4, N t = 6 and 
N t = 8. The temperature interval ranged from 4 • Tq upto 3 • 10 7 • Tq- Our results are shown on 
figure |[ A comparison can be done with the standard method, namely the result of [[Ij] for smaller 
temperatures. Here we present the results which are closest to the continuum limit, thus Nt = 
8. These results nicely follow the perturtbative predictions (see Jl0| , 11, [l2|], and also [13, 14]). 



However, more statistics is needed to determine the whole applicability region of the perturbative 
approach. It is important to point out, that at the highest temperature point, the pressure (within its 
statistical uncertainty) is already consistent with the Stefan-Boltzmann limit. 



7. Conclusion 

We presented two new methods to determine the equation of state for QCD. Compared to 
previous techniques, these methods allow to extend the temperature range by orders of magnitude. 
We developed a renormalization procedure of the pressure, for which no T = simulations are 
needed, and thus is less expensive than the methods, which were used previously. We also presented 
a new approach, which enables us to measure the pressure directly. Using this approach we gain 
a factor of N t in CPU time with respect to the standard method. Our approach can also be used to 
set the integration constant, which provides a reference point for the pressure. Based on these new 
methods we presented first results on the pressure of the pure gauge theory at temperatures that 
could not have been reached before. 
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Figure 5: The pressure, normalized to its Stefan-Boltzmann value, as a function of the temperature obtained 
by our new techique. Results with smaller discretization errors (N t = 8, blue circles), seem to fit improved 
perturbation theory, and also reproduce results obtained by the standard method at lower temperatures. At 
the highest temperature, 3 ■ 10 • Tq, the pressure (within its statistical uncertainty) is consistent with the 
Stefan-Boltzmann limit. 

[2] CP-PACS collaboration, M. Okamoto et al., Equation of state for pure SU(3) gauge theory with 
renormalization group improved action, Phys. Rev. D60 (1999) 094510 [hep- 1 at / 9905005]. 

[3] F. Karsch, E. Laermann and A. Peikert, The pressure in 2, 2+1 and 3 flavour QCD, Phys. Lett. B478 
(2000) 447-455 [hep-lat/0002003]. 

[4] MILC collaboration, C. Bernard et al., QCD equation of state with 2+1 flavors of improved staggered 
quarks, Phys. Rev. D75 (2007) 094505 [hep-lat/0611031]. 

[5] Y. Aoki, Z. Fodor, S.D. Katz and K.K. Szabo, The equation of state in lattice QCD: With physical 
quark masses towards the continuum limit, JHEP 01 (2006) 089 [hep-lat/ 0510084]. 

[6] F. Karsch, Recent lattice results on finite temperature and density QCD, In proceedings of XXV 
International Symposium on Lattice Field Theory (2007) PoS (LATTICE 2007) 015| . 

[7] J. Engels, J. Fingberg, F. Karsch, D. Miller and M. Weber, Nonperturbative thermodynamics ofSU(N) 
gauge theories, Phys. Lett. B252 (1990) 625-630. 

[8] S. Necco and R. Sommer, The N(f) = heavy quark potential from short to intermediate distances, 
Nucl. Phys. B622 (2002) 328-346 [hep-lat/0108008]. 

[9] M. Gockeler et al., A determination of the Lambda parameter from full lattice QCD, Phys. Rev. D73 
(2006)014513 [hep-ph/0502212]. 

[10] E. Braaten A. Nieto, Free Energy of QCD at High Temperature, Phys. Rev. D53 (1996) 3421-3437 

[hep-ph/9510408]. 

[11] K. Kajantie, M. Laine, K. Rummukainen and Y. Schroder, The pressure of hot QCD up to g**6 
ln(l/g), Phys. Rev. D67 (2003) 105008 [hep-ph/0211321]. 

[12] M. Laine and Y. Schroder, Quark mass thresholds in QCD thermodynamics, Phys. Rev. D73 (2006) 
085009 [hep-ph/ 60304 8]. 

[13] J. Blaizot, E. Iancu and A. Rebhan, Thermodynamics of the high-temperature quark gluon plasma 
(2003) [hep-ph/ 03031 85]. 

[14] J.O. Andersen, E. Braaten, E. Petitgirard and M. Strickland, HTL perturbation theory to two loops, 
Phys. Rev. D66 (2002) 085016 [hep-ph/0205085]. 



7 



